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ABSTRACT 



In this paper we consider the effects of resonance and electron scattering on the escape of Lyman a photons during cosmological 
hydrogen recombination. We pay particular attention to the influence of atomic recoil, Doppler boosting and Doppler broadening 
using a Fokker-Planck approximation of the redistribution function describing the scattering of photons on the Lyman a resonance 
of moving hydrogen atoms. We extend the computations of our recent paper on the influence of the 3d/3s-ls two-photon channels 
on the dynamics of hydrogen recombination, simultaneously including the full time-dependence of the problem, the thermodynamic 
corrections factor, leading to a frequency-dependent asymmetry between the emission and absorption profile, and the quantum- 
mechanical corrections related to the two-photon nature of the 3d/3s-ls emission and absorption process on the exact shape of the 
Lyman a emission profile. We show here that due to the redistribution of photons over frequency hydrogen recombination is sped up 

by AA'c/A'c 0.6% at z ~ 900. For the CMB temperature and polarization power spectra this results in |AC//C/| ~ 0.5% - 1% at 

/ > 1500, and therefore will be important for the analysis of future CMB data in the context of the Planck Surveyor, Spt and Act. The 

main contribution to this correction is coming from the atomic recoil effect (AN^/Nc 1.2% at z ~ 900), while Doppler boosting 

and Doppler broadening partially cancel this correction, again slowing hydrogen recombination down by AN^/N^ ~ 0.6% at z ~ 900. 
The influence of electron scattering close to the maximum of the Thomson visibility function at z ~ 1 100 can be neglected. We also 
give the cumulative results when in addition including the time-dependent correction, the thermodynamic factor and the correct shape 
of the emission profile. This amounts in AA'e /We ~ -1.8%atz~ 1160 and |AC;/C/| ~ 1% - 3% at / > 1500. 

Key words. Cosmic Microwave Background: cosmological recombination, temperature anisotropies, radiative transfer 



1. Introduction 



Motivated by 
Planck surveyor, 



the great experimental 
Spt and Act several 



prospects with the 
independent groups 



(e.g. iDubrovich & GrachevI 120051: IChluba & Sunvaevl 2006 : 



Kholupenko & Ivanchik 120061: iRubino-Martm et alj 12006 : 
Switzer & Hiratal 120081: Iwong & Scottll2007l) have investigated 



details in the physics of cosmological recombination and their 
impact on the theoretical predictions for the cosmic microwave 

■ background (CMB) temperature and polarization power spectra. 

' The declared goal for our theoretical understandin g of the 
ionization history is the ~ 0.1% accuracy level (e.g. see lHu et al] 
Ll995: Seljak et al. 200^ close to the maximum of the Thomson 
visibiht y function dSunvaey & Zeldovichlll970l) at z ~ 1100 
(e.g. see lSunvaev & Chluball2008l: iFendt et alj|2008l for a more 
detailed overview of the different previously neglected physical 
processes that are important at this level of accuracy). 

This paper is a continuation of our recent work on cosmo- 
logical recombination, in which we studied the effects of 3d- Is 
and 3s- 1 s two-photon processes on the dy namics of hydrogen re- 
combination (Chluba & Sunyaevll2009al) . Here we now wish to 
give the results for the changes in the Lyman a escape probabil- 
ity and free electron fraction when in addition accounting for the 
effects partial frequency redistribution related to resonance scat- 
tering of moving neutral atoms and electron scattering during 



this epoch. In our previous work we neglected this aspect of the 
problem, although in the standard textbook formulation based 
on a Fok ker-Planck exp ansion of the frequency redistribution 
function ( lRybickill2006ll we obtained these results already some 
time ago. Here we explain the main r esults of these computations 
which we also partly used elsewhere jRubino-Martfn et al .120081: 
IChluba & Sunyaevll2008l) . and also refine our computations in- 
cluding the 3d- Is and 3s- Is two-photon corrections. 



Send offprint requests to: J. Chluba, 
e-mail: jchluba@mpa-garchir\g . mpg . de 



It is well known (e.g. see lRvbicki & deirAntoniolI994l) that 
for the conditions in our Universe (practically no collisions) the 
frequency redistribution function for photons scattering off mov- 
ing ato ms is given by th e so called type-II redistribution as de- 
fined in lHummeii ( Il962h . The main physical processes which are 
accounted for in the Fokker-Planck expansion of this frequency 
redistribution function are due to (i) atomic recoil, (ii) Doppler 
boosting, and (iii) Doppler broadening. All three physical pro- 
cesses are also well-known in connection with the Kompaneets 
equation which describes the repeated scattering of photons by 
free electrons. Atomic recoil leads to a systematic drift of pho- 
tons towards lower frequencies after each resonance scattering. 
This allows some additional photons to escape from the Lyman 
a resonance and thereby speeds hydrogen recombination up, a s 
already demonstrated earlier by iGrachev & DubrovichI (|2008|) . 
We found very similar results for this process some time ago 
(e.g. see footnote 10 in Chluba & Su nvaev..2009bi) . which here 
we shall present in detail and also refine including additional 
corrections simultaneously. 
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However, in the analysis of lGrachev & Dubrovichl ( l2008l) the 
effect due to (ii) and (iii) were not taken into account. Like 
atomic recoil Doppler boosting leads to a systematic motion 
of photons, but this time towards higher frequencies. Therefore 
it is expected to slow recombination down. In contrast to this 
Doppler broadening can lead to both an increase or a decrease 
in the escape probability depending on where the photon ini- 
tially is emitted. As we explain here, if the photons are initially 
emitted in the vicinity of the Doppler core line diffusion helps 
to bring some of them towards the red wing, before they actu- 
ally die (mainly due to two-photon absorption to the third shell). 
Similarly, for photons emitted on the blue side of the resonance 
line broadening allows some finite number of them to transverse 
the Doppler core. In the no line scattering approximatiorQ this 
would not be possible, so that in both case the escape fraction is 
increased. In contrast to this, for photons emitted on the red side 
of the resonance the effect of Doppler broadening decreases the 
escape fraction, since even up to ~ 100 Doppler width below the 
line center a significant fraction of the photons still returns close 
to the Doppler core, where they die efficiently. As we show here, 
the combination of Doppler boosting and Doppler broadening in 
total leads to an additional decrease in the escape probability as 
compared to the no line scattering approximation. 

For the expected correction due to electron scattering very 
similar arguments apply. However, there are some important dif- 
ferences: (i) electron scattering is expected to become less im- 
portant at lower redshifts, since the free electron fraction de- 
creases with time; (ii) in contrast to resonance scattering for 
Lyman a photons the electron scattering cross section is achro- 
matic; and (iii) due to the smaller mass of the electron the recoil 
effect is ~ 2000 larger Nevertheless, it turns out that during hy- 
drogen recombination electron scattering can be neglected in the 
analysis of future CMB data. This is because of its much smaller 
cross section in comparison with line scattering and the decreas- 
ing number density of free electrons (see Sect. 12.21 ). 

We would like to mention that while this paper was in prepa- 
ration another investig ation of this problem was carried out by 
iHirata&ForbesI (120091) . The results obtained in their work seem 
to be in good agreement with those presented here. 

2. Additions to tlie kinetic equation for tlie photons 
in the vicinity of the Lyman a resonance 

Here we give the additional terms for the photon radiative 
transfer equation which are necessary to describe the effect 
of resonance and electron scattering in the Lyman a escape 
problem during cosmological hydrogen recombina t ion. We will 
use the same nota t ion as in IChluba & SunvaevI (l2009bl) and 
IChluba & SunvaevI ( l2009ah . also introducing the dimensionless 
frequency variable x - v/(l + z) and photon distribution, .^V., - 
Nx/(i - Ny/(l +z)^, with Ny - ly/hv, where ly is the phys- 
ical specific intensity of the ambient radiation field. The photon 
occupation number then is ny - lyjlhv^ . Note that with this 
choice of variables the redshifting of photons due to the Hubble 
expansi on is automatically taken in to account in x (for more de- 
tails see lChluba & Sunvaevll2009bh . 

It is clear that Lyman-cy line and electron scattering (both in- 
cluding the Doppler-broadening, recoil and induced scatterings) 
only lead to the redistribution of photons over frequency, but do 
not change the total number of photons in each event. Also a 



' In this approximation only true line emission and line absorption 
and redshifting of photons are included in the computation. The redis- 
tribution of photons over frequency is neglected. 



blackbody spectrum with Ty - should not be altered by these 
processes. Within the Fokker-Planck formulation of the corre- 
sponding processes these requirements are directly fulfilled. 

2.1. Lyman-a resonance scattering 

The contribution to the collision term due to redistribution of 
photon by resonance scatt ering off moving atoms can be written 
as (e.g. see lRvbickill2006l) 

C[Ny-\ I, = J 5R(v, V)Ny, (1 + «,) dv' 

- ^ niv', v)Ny(\ + ny.) dv' , (1) 

where 'R{v, v') is the frequency redistribution function for the 
scattering atom, which for conditions in the Universe (practically 
no electron or proton collisions!) is pur ely due to the Do ppler ef- 
fect (type-II redistribution as defined in ^HummeJ I 1 9621) . 

As shown in iRvbicki ( 2006i) . within a Fokker-Planck formu- 
lation for the case of Doppler redistribution Eq. ([TJ can be cast 
into the form 



C[Ny]\, ^ p'j, (r,Ni, 



ov 



d_Ny_ _h_N^l, c^N^ 
dv v2 ^ kT, v2 ^ 2v2 



(2) 



where cr.. 



^ ^ denotes the resonant scattering cross 
section and Avd the Doppler width of the Lyman a reso- 
nance. The first term in brackets (oc dyny) describes the com- 
bined effect of Doppler boosting (it is of the order ~ V~/c~, 
where V is the velocity of the atom) and Doppler broadening 
(~ V/c), while the second term (oc ny[l + «,,]) accounts for 
atomic recoil (~ /zv/whc^) and stimulated scatterings. Following 
iRvbicki & deirAntoniol ( ll994l) we have used the diffusion coef- 
ficient D oc (pviv), where (pv{v) is the normal Voigt-profile. We 
wil l neglect c orrections due to non-resonant contributions (e.g. 
see iLeel2005|) in the scattering cross section, which would lead 
to a different frequency dependence far away fro m the resonanc e 
(e.g. Rayleigh scattering in the distant red wing ( |Jacksonll998l) ). 

It is important to note that Eq. (|2]i simultaneously includes 
the effects of line diffusion, atomic reco/Q and stimulated scat- 
tering In this formulation it therefore preserves a Planckian 
photon distribution Ny^ with Ty s T^.. This can be easily 

verified when realizing that 4-^ s -tt^-VII + -r^l Also 

® Oy V- kly y- \ 2v- J 

one can easily verify that in the Fokker-Planck formulation 
C[A^v]lr dv = 0. 
In Eq. ^ we also took into account the fact that not ev- 
ery scattering leads to the reappearance of the photon, since per 
scattering the fraction 1 -pic of photons disappear in other chan- 
nels, i.e to higher levels and the continuum. Here pj is the sin- 
gle scattering albedo which in our formulation is equivalent to 



the one photon emission probability p'^. However, since p'J i s 
always very close to unity (e.g. see IChluba & Sunvaevll2009a) . 
one could also neglect this detail here. 



This terms was first introduced bv lBaskol jl978l.ll981h 
^ For the escape of Lyman a photons during hydrogen recombination 
this term is not important. 
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Fig. 1. Time evolution of the photon distribution for single (narrow-line) injection at the line center. The death probability for a 
3-shell hydrogen atom was used and electron scattering has been neglected. We use the time-variable t = J c cr^N\{ At. 



The corresponding term in the variables x and A^^ then reads 



Irj- 



dx 



X (pyiy) 



dx x^ X- \ 2x^ 



(3) 



where we have made the substitutions rj = (1+z), ^ - ^ and 

v = 41 Note that = x?j where X2i = V2i/[l+z] and 
mil is the mass of the hydrogen atom. This term has to be added 
to the radiative transfer equation which inclu des the effect of line 
emission and absorption and can be found in lChluba & SunvaevI 
(l2009b|) for the normal ' 1 + 1 ' phot on formulation of the problem 
and in Tchluba & SunvaevI (l2009al) for the two-photon formula- 
tion. 



2.2. Electron scattering 

The contribution to the collision term due to scattering off free, 
non-relativistic electrons can be described with the Kom paneets- 
equat ion. Due to the similarity with Eq. Q (see also iRvbickil 
l2006l) it is straightforward to obtain the corresponding terms for 
our set of variables: 



C[N,]\^^(nN,e, — x'' 



d N, N, 
dx x^ x^ 



1 -H 



2x^ 



(4) 



where crj ~ 6.65 x 10""^ cm~ is the Thomson cross section and 
0e - k^T^jm^c^. Again one can clearly see that the electron scat- 
tering term preserves a Planckian photon spectrum for Ty = Tg, 
and that / C[A^v]lr dv = 0. 

2.2.1. Relative importance of electron scattering 

Since the line-profile (py is a strong function of frequency, res- 
onance scattering is most important close to the Lyman a line 
center, while in the very distant wings electron scattering is ex- 
pected to dominate. Comparing the diffusion coefficients in fre- 



quency space for resonant and electron scattering 



a.9x 10" 



Nu 0v(v) 



(5a) 



~ 2.1x10-^1 +z)-fx2, (5b) 

^Vls 

shows that at redshift z ~ 1400 (where Ne/N\s ~ 4) in the line 
center resonance scattering is ~ 2.0 x 10^ times more important 
than electron scattering, and only at \x£,\ > 100 Doppler width 
electron scattering is able to compete with line scattering. 

Due to the changes in Ne/Nu the ratio (|5]l is a strong func- 
tion of redshift. However, electron scattering is expected to in- 
fluence the evolution of photons close to the line center signifi- 
cantly only at redshifts z > 2500, i.e. well before the main epoch 
of hydrogen recombination. Therefore one expects that electron 
scattering has a small impact on the development of the photons 
close to the center of the Lyman-a transition and hence on the 
escape probability during hydrogen recombination. 

3. Illustrative time-dependent solutions for different 
initial photon distributions 

In this Sect, we illustrate the main physical effects related to res- 
onance scattering and electron scattering. For this we numer- 
ically solved the radiative transfer equation injecting a single 
narrow-line at different distances from the line center. For the 
computations we include the frequency redistribution of pho- 
tons, redshifting a nd real ab sorption using the normal ' 1 + 1 ' pho- 
ton picture (see iChluba & Sun vaev 2009b). We neglect the ef- 
fects due to two-photon corrections here. Furthermore, we shall 
assume that the solution for the electron number density and 
the 1 s-population a re given by the output of the Recfast code 
(ISeager et al.lll999l) . A few words about the PDE-solver can be 
found in the Appendix lAl 

3.1. Time-dependent solutions 

In Fig. [T]we present the results for single injection of photons 
at the Lyman-a line center In practice we use a Gaussian ini- 
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. 2. Time evolution of the photon distribution for single (narrow-line) injection on the red side of the Lyman-a resonance at 
different distance from the line center. The death probability for a 3-shell hydrogen atom was used and electron scattering has been 
glected. We use the time-variable t - J c crrNu df . 



Chluba and Sunyaev: Cosmological Hydrogen Recombination: influence of resonance and electron scattering 5 




Fig. 3. Time evolution of the photon distribution for single (narrow-line) injection at xd = +10. The death probability for a 3-shell 
hydrogen atom was used and electron scattering has been neglected. We use the time-variable t - j c o-^Nu dt. 



tial photon distribution which is centered at the injection fre- 
quency XD.i and has a width cr^ ~ 5 x 10"~. Furthermore, we 
re-normalized by a convenient factor such that induced effects 
are negligible. We started our computation at injection redshift 
Zi - 1400, i.e. close to the time where the maximum of the 
CM B spectral distortion due to the Lyman-a transition appears 
jRubifio-Martfn et al. 2006). At this redshift roughly 20% of all 
hydrogen atoms have already recombined and the death proba- 
bi hty for a 3-shell hydrogen atonfl is pd ~ 5.6 x 10""* (see Fig. 1 
in lChluba & Sunvaev||2009bh . 

From Fig.[T]one can see that after a short time the initial pho- 
ton distribution has broadened significantly, bringing photons to 
the wings of the Lyman-a transition. After t - j ccr^ Nu dt ~ 
10"* the death of photons in the line center becomes important, 
owing to the fact that ~ 10""* is so small. The solution remains 
very symmetric until r ~ few x 10^° and only then redshifting 
due to the expansion of the Universe starts to become important 
(as we will see line-recoil only affects the photon distribution at 
the level of few percent in addition). When the bulk of photons 
reaches a distance xu 100 still a sizable amount of them re- 
mains on the blue side of the Lyman-a line, and only when the 
maximum of the photon distribution reaches xo ~ - 150 the evo- 
lution starts to become dominated by redshifting and absorption 
only, with very small changes because of frequency redistribu- 
tion. 

In Fig. |2] we present the results for single injection of pho- 
tons at different distances to the line center. Again photons were 
injected at Zi = 1400. Focusing on the case XD,i = -10, one can 
again observe the fast broadening of the initial photon distribu- 
tion. However, now the characteristic time for line scattering has 
increased by a factor of ~ 2.3 x 10^ because frequency redis- 
tribution akeady takes place in the wings of the Voigt-profile. It 
is important to note that due to line scattering photons strongly 
diffuse back into the line center and thereby increase the possi- 
bility of being absorbed. Also one can see that due to diffusion 



The main contribution to the death of photons is due to the two- 
photon absorption to the 3d-state. Including more sheUs the death 
probabihty changes by les s that 10% during hydrogen recombination. 
ichluba & Sunvaevll2009 j) . 



some photons even reach far into the blue side of the Lyman-a 
resonance. Again only after the bulk of photons has reached a 
distance of ~ -150 redshifting and absorption play the most 
important role for the evolution of the photon distribution. 

Looking at the other two cases, it becomes clear that for in- 
jection at xo.i = -50 still a few photons do diffuse back to the 
line center, whereas for xo.i - -100, practically all photons re- 
main below xu ~ -50 at all times. Comparing the maxima of 
the final photon distribution (at t ~ 2.5 x 10") for all the dis- 
cussed cases shows that as expected the efficiency of absorption 
decreases when increasing XD,i. 

It is also interesting to look at cases when injecting photons 
on the blue side of the Ly man-resonance. In this case all pho- 
tons have to pass at least once through the resonance before they 
can escape and one expects that many photons die during this 
passage. In Fig. |3] we show the results for single injection at 
xd = +10. At the beginning the evolution of the spectrum looks 
very similar (except for mirror-inversion) to the case of injection 
at xd = -10. However, at late times one can see that the amount 
of photons reaching the red side of the Lyman-a resonance is 
significantly smaller Indeed this amount is comparable to the 
case of injection directly at the center. 

3.2. Escape probability for single narrow line injection 

Given an initial photon distribution one can compute the total 
number of photons that survive the evolution over a period of 
time for the given transfer problem. Here we assume that only 
at time t - fresh photons are appearing. Comparing the total 
number of photons at the final stage with the initial number then 
yields the numerical escape or survival probability for the given 
diffusion problem 



Pesc(Zi,Zf) 



Ny(Zf) 

NJzi) 



1+Zi 



1 +Zf 



/ NAzf)dx 
jNAzi)dx' 



(6) 



where Ny{z) - J Nyiz) dv is the number density of photons at 
redshift z. The factors (1 + z)^ account for the changes in the 
scale factor of the Universe between the initial and final redshift. 
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Due to the expansion of the Universe photons redshift to- 
wards lower frequencies. Neglecting any redistribution process, 
with time this will increase the distance of the initial photon dis- 
tribution to the line center and thereby decrease the probability 
of real line absorption. Assuming that the initial photon distribu- 
tion is given by a ^-function then with Eq. ^ one obtains 

for this case. Here Tabs is the absorption optical depth between 
the initial redshift Zi and z- 

We now want to compare the differential escape probability 
Eq. (|2l) with the numerical results obtained when also including 
the redistribution of photons over frequency. The results of the 
the previous Section suggest the following: 

(i) For photons injected close to the line center the diffusion 
due to resonance scattering helps to bring photons towards 
the wings. In comparison to the case with no scattering this 
should increase the escape probability. 

(ii) At intermediate distances on the red side of the line center 
(jcd ~ -50 to -100 Doppler width) line diffusion brings 
some photons back to the Doppler core and thereby should 
decrease the escape probability in comparison to the case 
without line scattering. 

(iii) Far in the red wing of the line (xo ^ -100) the escape frac- 
tion will depend mainly on the death probability and the ex- 
pansion rate of the Universe. In this regime line scattering 
does lead to some line broadening, but should not affect the 
escape probability significantly anymore. 

(iv) The escape probability for injections on the blue side of the 
resonance becomes nearly independent of the initial distance 
to the line center and should be comparable to the one inside 
the Doppler core. 

It is easy to check these statements numerically. For this we 
performed a sequence of computations injecting photons at dif- 
ferent distances from the line center and following their evo- 
lution until the initial maximum of the photon distribution has 
reached XD,t - XD,i - xu,s- We then computed the escape or sur- 
vival probability as defined by Eq. ^ for the given diffusion 
problem as a function of the injection frequency, xo.i, injection 
redshift, Zi, and termination redshift, Zt, which directly depend^ 
on the value of xu,s- 

Since the absorption cross section in wings of the line scales 
like oc I/jcd, even beyond xo ~ -10^ still percent-level absorp- 
tion can occur, which should be taken into account when com- 
puting the total escape probability until redshift Zt - 0. However, 
the effect of resonance scattering becomes negligible at this dis- 
tance from the line center (see below) and the time evolution 
in principle can be described fully analytically. For simpUcity 
we neglected this additional complication and typically chose 
•^D,s ~ 10"*, which ensured that the remaining absorption will 
only lead to modifications of AP/P < 10^^ to the obtained es- 
cape probability. Up to this level of accuracy, the obtained curves 
presented in this Section can be considered as the frequency- 
dependent total escape probability until Zt = 0. 

In Fig. ID we present some results for computations of the 
frequency-dependent escape probability, Pesc{xD,Zi,Zt), for in- 
jection redshifts Zi = 1100 and 1300. For comparison we also 
give the corresponding escape probabilities, Pf^c"' - e^^'^', 
Eq. (|7l), for ^-function injection when neglecting line scattering. 
At large distance (xd < -150) from the line center Pesc practi- 
cally coincides with P^^s^^ in all presented cases. As mentioned 

' For simplicity we used Zi = + ^D,tAyD(Zi)/v2i]. 




Fig. 4. Escape probability, Pef,c(xu, Zi,Zi), for single (narrow line) 
injection at different distances from the line center and initial 
redshifts z\ - 1100 and 1300. The death probability for a 3- 
shell hydrogen atom was used and electron scattering has been 
neglected. For the given curves we set xd,s - 10"*, such that 
Pe!ic{xD,Zi,Zi) ~ PescixD,Zi,Q)- For Comparison also the analytic 
result, Pf'sc'N for 5-function injection including only pure absorp- 
tion without line scattering is shown. 

above this behavior is expected since line scattering should not 
strongly affect the evolution of the line anymore. At intermediate 
distances from the line center the inclusion of line scattering in- 
deed decreases the escape probability in comparison to the cases 
without scattering. Looking in detail at the dependence of Pesc 
close to the center of the line shows that the presumptions (i) 
and (iv) also hold. Our computations clearly show that there is a 
non- vanishing escape probability for photons from the blue side 
of the line, which in the case of pure absorption is practically 
zerc0. This probability is nearly constant extending even into the 
core of the line and down to xq 2. 

* There is a small difference close to the line center due to the fact 
that we used i5-function injection for the computation of Pefc"* instead 
of the Gaussian that was used in the numerical computation. However, 
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Fig. 5. Frequency xd < at which the modification due to line 

scattering becomes e percent. Here AP/P = [PgL ^^esc'^]/^e's?'^ 
where for P^^ line scattering was taken into account. The death 
probability for a 3-shell hydrogen atom was used and we set 
.^D s - lO'* for all curves. 



In order to understand up to which distance to the line cen- 
ter the effect of resonance scattering is important we compared 
the results for the escape probability including line scattering 
with the analytic no-scattering solution, asking the question at 
which distance in the red wing the modification due to line scat- 
tering becomes e percent. In Fig. |5] we summarize the results 
of such comparison. It is clear that at all redshifts of interest 
line-scattering is only important for xo > -few x 10^, but at 
the percent-level in principle may be neglected below this fre- 
quency. We made use of th is result already in some earlier works 
(IChluba & Sunvaevll2008l) . 

3.2.1. Role of atomic-recoil 

Every resonance scattering due to atomic-recoil leads to a 
small shift of the photon energy towards lower frequencies. 
The strength of the recoil due to the frequency-dependence of 
the scattering cross section is a strong function of photon en- 
ergy, peaking close to the Lyman-a line center, and dropping 
rather strongly in the damping wings. This is in stark contrast to 
electron-recoil, for which the scattering cross section is practi- 
cally independent of frequency. 

To understand the importance of the atomic recoil effect for 
the differential escape probability we therefore performed sev- 
eral computations of the frequency-dependent escape probability 
for injection of photons at different distances from the line cen- 
ter explicitly neglecting the effect of atomic recoil. In Fig.|6]we 
present the correction to the escape probability which is only due 
to the atomic-recoil term. As expected, atomic recoil helps pho- 
tons to escape in the whole range of frequencies. However, due 
to the decrease in the scattering cross-section the corresponding 
correction becomes very small at distances below xu ~ -100 to 
-150. Also the amplitude of the effect increases towards lower 
redshifts, simply because more hydrogen atoms have become 
neutral. The largest correction is coming from the line center 



this will only make the transition P^fc ' 
will not change the main conclusion. 



4 

a, 3 



= 1100 
= 1200 
= 1300 
= 1400 
= 1500 



0L_ 
-150 



-100 



-50 



Fig. 6. Relative difference in the escape probability for single 
(narrow line) injection at different distances from the line center 
when including the effect of atomic recoil. Here AP/P = [Pesc - 
■Pesc"'"]/-Pesc"'"' whcrc for P™"'"'"' the term due to atomic recoil 
was neglected. The death probability for a 3-shell hydrogen atom 
was used and we set xus - 10"*. 



a. 




Fig. 7. Relative difference in the escape probability for single 
(narrow line) injection at different distances from the line cen- 
ter when including electron scattering. Here AP/P = [Pl^^ - 
^'escl/^'esc, whcrc for Pgsj, electron scattering was taken into ac- 
count. The death probability for a 3-shell hydrogen atom was 
used and we set xd,s = 10"* for all curves. 



less steep, but otherwise 



and in practically constant over the whole Doppler core and the 
blue side of the resonance. We will see below that the total cor- 
rection to the Lyman a escape probability is very similar to the 
value obtained for injections close to the line center (see Sect.|4]i. 
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3.2.2. Role of electron scattering 

In Fig.|7]we show the relative difference in the escape probabil- 
ity for single (narrow line) injection at different distances from 
the line center when including electron scattering. As expected 
electron scattering has a similar effect as resonance scattering, 
helping photons to escape more efficiently from the line cen- 
ter, but bringing some photons from the wings back into the 
Doppler-core, diminishing the probability of their survival. At 
higher injection redshift the differences become larger, due to 
the increase of the number of free electrons as compared to the 
number of neutral hydrogen atoms. At Zi < 1200 the relative dif- 
ference becomes smaller than ~ 1 .2% in the whole considered 
range of injection frequencies. Close to the maximum of the vis- 
ibility function Zi ~ 1 100 one does not expect a large correction 
due to electron scattering. In addition it is clear that the increase 
in the escape in the Doppler-core should be partially canceled 
by the decrease in the red wing. As we will see in Sect.|4]the net 
effect of electron scattering on the Lyman a escape probability 
during hydrogen recombination is always < 1% at z < 1600. 



4. Changes in the Lyman a escape probability 
during hydrogen recombination 

In this Section we now present the results for the changes in 
the Lyman a escape probability during hydrogen recombination. 
Our approach here is very similar to the one used in our ear- 
lier, semi-analytical works dChluba & Sunvaevll2009b ia!). Given 
the solution for the populations of the different hydrogen lev- 
els we numerically solve the transfer equation for the Lyman 
a problem obtaining the spectral distortion in the vicinity of the 
Lyman a resonance at different redshifts. From this we can com- 
pute the effective escape probability by convolving this distor- 
tion with the corresponding Lyman a absorption profile. We also 
followed a very similar approach in our previous computations 
of the radiative transfer problem during helium recombination, 
w here some of the resuhs obta ined in that case were already used 
in lRubino-Martm eta n (l2008h . 

We will start by discussing the results in the standard ' 1 + 1 ' 
formulation (Sect. 14.11 ). We then include the effect due to the 
thermodynamic correction factor fy (Sect. 14.2b , which was in- 
troduced earlier ( IChluba & Sunvaevl2009bllai) using the detailed 
balance argument. Finally we shall also include the corrections 
to the 3d- Is and 3s- Is two-photon emission profile (Sect. |43] l. 

4.1. Results in the standard 'I + I' photon formulation 

In Figure[8]we present the results for the escape probability using 
the standard '1 -H 1' photon formulation. In this case the emis- 
sion and absorption profile are given by the normal Voigt profile. 
We also included the full time-dependence of the problem in the 
computations of the line emission rate and th e absorption optical 
depth. In the no scattering approximation (IChluba & SunvaevI 
I2009bh this leads to the dashed curve shown in Fig. [8] 

As mentioned earlier (IChluba & SunvaevI l2009bh the stan- 
dard '1 + r photon formulation has several discrepancies, i.e. 
leading to an unphysical self-feedback of Lyman a photons at 
low redshifts {z < 800 - 900). Nevertheless, one can study the 
influence of the redistribution of photons by resonance and elec- 
tron scattering even in this approach and as we will see one ob- 
tains very similar results for the effect of resonance scattering in 
comparison with the more complete formulation using the two- 
photon picture (Sect. |43] ). 
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Fig. 8. Changes in the Lyman a escape probability for the stan- 
dard ' 1 -I- 1 ' photon formulation. The dashed l ine shows the result 
obtain ed in the no scattering approximation dChluba & SunvaevI 
I2009ah . 
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Fig. 9. Changes in the Lyman a escape probability due 
to the thermodynamic correction factor The dashed line 
shows the result obtained in the no scattering approximation 
(IChluba & Sunvaevl2009al) . 



In Figure [8] we show the separate correction due to atomic 
recoil (thin solid line). We obtained this curve by taking the dif- 
ference of the escape probabilities for the case with all correc- 
tions due to line and electron scattering included and the one in 
which line recoil was switched off. The importance of recoil in- 
creases towards lower redshifts reaching the level of AP/P ~ 6% 
at z ~ 600. If we look at the results presented in Fig.|6|for the 
case of single narrow line injection, we can even see that the 
total recoil correction seen in Fig.|8]is very close to the value ob- 
tained for line center injection. This is expected, since the largest 
contribution to the total value of the escape probability always 
comes from the Doppler core. 
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We can also see that the effect of electron scattering (dotted 
curve) is very small, leading to a correction < 1% at z < 1800. 
Close to the maximum of the Thomson visibility function the ef- 
fect of electron scattering is negligible. This curve was computed 
using the numerical results in which we switched off electron 
scattering and then compared it to the one where it was included. 

Finally, we also computed the contribution that can be at- 
tributed to the effect of Doppler boosting and Doppler broaden- 
ing (dash-dotted curve). For this we computed the escape proba- 
bility when neglecting electron scattering and atomic recoil, but 
only including the line diffusion term. We then took the differ- 
ence to result obtained in the no scatter ing approximation, as 
given earlier dChluba & Sunvaevll2009bh . One can see that the 
diffusion term results in a decrease of the escape probability at 
low redshifts. However, this decrease is about 3 times smaller 
than the increase in the escape probability due to atomic recoil. 
Therefore the net effect due to resonance scattering is an increase 
in the escape probability, reaching AP/f ~ -1-4% at z ~ 600. As 
explained in Sect. 13.21 this shows that the decrease in the red 
wing escape probability due to the return of photons towards 
the Doppler core by line diffusion is more important than the 
increase of the escape fraction from within the Doppler core 
caused by Doppler broadening. 

We would like to mention that the small variability in the dif- 
fusion contribution at z ~ 600 is likely due to some details in our 
numerical treatment. However, we expect that the corresponding 
result is converged at the ~ 10% level of the correction, which is 
sufficient for our purposes here. 

4.2. Effect of the thermodynamic corrections factor 

If we now in addition include the frequency-dependent asym- 
metry between the emission and absorption profile due to the 
therm odynamic correction factor which was introduced earlier 
(I Chl uba & Sunvaev 2009b a), we obtain the results presented in 
Fig- HI The dashed lin e again shows the correctio n in the no scat- 
termg approximation (IChluba & Sunvaevl2009al) . The main cor- 
rection do to the redistribution of photons over frequency again 
is due to the line recoil term (thin solid line). One can see that it 
is practically the same as in the previous case (see Fig.[8]l. Also 
the total correction due to electron scattering did not change very 
much. In both cases the difference was smaller than ~ 5% on the 
correction. However, the correction due to the line diffusion term 
seems to be slightly increased, suggesting a fy induced correc- 
tion to the correction that is not completely negligible. 

4.3. Corrections due to the shape of the emission profile 

Finally, we also ran the code including th e correct shape of 
the 3d/ 3s- Is emission and absorption profile dChluba & SunvaevI 
l2009ah . The results of these computations are shown in Fig.fTOl 
The dashed line ag ain shows the correction i n the no scatter- 
ing approximation (IChluba & Sunvaevll2009al) . The dotted line 
in addition indicates the correction that was associated with the 
effec t of the emission profile in the no scattering approximation 
dChl uba & Sunvaev 2009a). We also computed the pure profile 
correction using the numerical results obtained when including 
the redistribution of photons and obtained the dash-dotted curve. 
As one can see the difference to the no redistribution case is very 
small. Therefore we did not compute the pure recoil correction, 
the line diffusion correction or the correction due to electron 
scattering, since they should also be very similar to the contribu- 
tions shown in Fig.|9l 
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Fig. 10. Changes in the Lyman a escape probability due to the 
shape of the emission profile. The dashed li ne shows the result 
obtain ed in the no scattering approximation dChluba & SunvaevI 
I2009ah . 
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Fig. 11. Changes in the Lyman a escape probability due to 
the different processes under discussion here. The dashed line 
shows the result obtained in the no scattering approximation 
(IChluba & Sunvaevl2009al) . 



5. Corrections to the ionization liistory 

In this Section we now give the expected correction to the ioniza- 
tion history when including the processes discussed in this pa per 
For this we modified the Recfast code dSeager et al.lll999h . so 
that we can load the pre-computed change in the Sobolev escape 
probability studied here. 

In Fig. [TT]we present the final curves for AP/P as obtained 
for the different processes discussed in this paper In Fig. [12] 
we show the corresponding correction in the free electron frac- 
tion computed with the modified version of Recfast. The atomic 
recoil effect alone (thin solid line) leads to ANe/Ng ~ -L2% 
at z ~ 900. This is in very good agreement with the result 
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Fig. 12. Changes in the free electron fraction due to the 
different processes under discussion here. The dashed line 
shows the result obtained in the no scattering approximation 
(IChluba & Sunyaevll2009ah . 

of iGrachev & Dubrovic h' ('2008') . We already quoted this result 
earlier (see footnote 10 in Chluba & Sunyaevll2009bl) . however 
there we just estimated the change in the free electron frac- 
tion using our full numerical result for the recoil correction on 
the Lyman a escape probability, without running it trough the 
Recfast code. Including electron scattering and all terms (line 
recoil and the diffusion term) for the redistribution of photons by 
the Lyman a resonance we obtain the dotted line. Here the to- 
tal correction due to redistribution of photons now only reaches 
AA^e/A^e ~ -0.6% at z ~ 900. As we have seen in Sect, g] 
this is due to the fact that the diffusion term slow recombina- 
tion down again, since photons from the red wing return close 
to the Doppler core, where they die efficiently again. Finally, 
the total correction including all the effects of photon redis- 
tribution and the correction due to the time-dependence, ther- 
modynamic fa ctor and shape of the pro file, which were dis- 
cussed earlier dChluba & Sunvaevll2009ah . has a maximum of 
ANe/Ne ~ -1.8% at z ~ 1160. Here the main contribution is 
coming from the the time-dependent correction and thermody- 
namic factor as explained in Chluba & Sunyaev (2 009al) . 

In Fig.[T3]we finally show the changes in the CMB temper- 
ature and polarization power spectra. The corrections to ANe/Ne 
related to the redistribution of photons over frequency alone (up- 
per panel) results in changes to the TT and EE power spectra, 
with peak to peak amplitude ~ 0.5 % - 1 % at / > 1500. When also 
including the processes discussed in lChluba & SunvaevI (l2009al) 
at / > 1500 we find a cumulative correction of |AC//C/| ~ 1% 
for the TT power spectrum and |AC//C/| ~ 2% - 3% for the EE 
power spectrum. It will be important to take these changes into 
account in the analysis of future CMB data. 

6. Conclusions 

In this paper we have considered the effect of frequency redis- 
tribution on the escape of Lyman a photons during hydrogen 
recombination. We have shown that line recoil speeds hydrogen 

recombination up by ANe/Ne 1.2% at z ~ 900. On the other 

hand, the combined effect of Doppler boosting and Doppler 
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Fig. 13. Changes in the CMB temperature and polarization 
power spectra. The upper panel shows the changes due to the re- 
distribution of photons by line and electron scattering only. The 
lower panel shows the cumulative result in addition including 
the time-dependent correction, the thermodynamic factor, and 
the correction due to the shape of the em ission profile, as dis- 
cussed earlier jChluba & Sunvaevll2009al) . 



broadening at different distances from the line center slows hy- 
drogen recombination down by ANg/Ng ~ 4-0.6% at z ~ 900. 
As explained in Sect. |3] line diffusion (including both Doppler 
boosting and Doppler broadening) increases the escape fraction 
for photons that are emitted in the vicinity of the Doppler core in 
comparison with the value obtained in the no scattering approx- 
imation. In particular some small fraction of photons that are 
emitted on the blue side of the resonance can still escape, since 
due to line diffusion they pass through the Doppler core faster 
than dying there. On the other hand, for photons that are emitted 
at -fewx 10^ < xo < -10 (i.e. in the red wing) it becomes harder 
to escape, since line diffusion brings some of these photons back 
close to the Doppler core, where they are absorbed efficiently. 
For photons that are emitted at xq < -few x 10^ the redistri- 
bution over frequency can be neglected. We also showed that 
electron scattering has a minor effect on recombination dynam- 
ics at redshifts z < 1400. In total the redistribution of photons 
over frequency leads to a speed up of hydrogen recombination 
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by AA^e /Ne 0.6% at Z ~ 900 (cf. Fig. [T2|. This results in im- Sunyaev. R. A. & Zeldovich, Y. B. 1970, Astrophysics and Space Science, 7, 3 

portant changes to the CMB temperature and polarization power Switzer, E. R. & Hirata, C. M. 2008, Phys. Rev. D, 77, 083006 
spectra (see Fig. [l3]for details), which should be taken into ac- & Scott, D. 2007, MNRAS, 375, 1441 

count for the analysis of future CMB data. 

In addition, we would Uke to mention that the cu- 
mulative changes (inclu ding the processes discussed in 
IChluba & Sunyaev ' (2009a') and those of this work) in the Lyman 
a photon escape probability will be very important for precise 
co mputations of the cosmolog ical recombination spectrum (e.g. 
see ISunvaev & Chluball2007l for review and references). Here 
it is interesting that the changes in the shape of the recombina- 
tion lines connected with electrons passing through the Lyman a 
channel are expected to be ~ 10% at z ~ 1400 (in comparison to 
~ 2% for Ne at z ~ 1200). Observing the cosmological recom- 
bination lines and looking at their exact shape would therefore 
provide a more direct and -4-5 times more sensitive probe for 
the physics of cosmological recombination than with the CMB 
temperature anisotropics. 

Appendix A: Computational details 

A.1. Solver for the differential equations 

In order to solve the photon transfer equation we used the solver 
D03PPF from the NAcQ-Library. It provides possibilities for ex- 
tensive error control and adaptive remeshing. In particular for 
computations with narrow initial spectra or low line scattering 
efficiency this feature became very important. However, remesh- 
ing also leads to an additional loss of accuracy for long integra- 
tions and therefore has to be applied with caution. 

Typically we used ~ 2500 - 5000 grid-points for the repre- 
sentation of the photon distribution and required relative accura- 
cies e ~ lO"*" - 10"^. We checked the convergence of the results 
by varying the accuracy requirements and number of grid-points, 
and also by running several test problems for which analytic so- 
lutions exist. 
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